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LONG-TERM  GOALS 

The  long-term  objective  is  to  develop  knowledge  of  factors  controlling  the  mean  and  the  variability  of 
acoustic  fields  in  the  shallow  ocean  environment,  in  the  frequency  band  50  Hz  to  3000  Hz.  Propagation 
variability  is  an  inescapable  complicating  factor  for  both  active  and  passive  sonar  systems,  and  for 
underwater  acoustic  communications.  Reliable  predictions  of  temporal  and  spatial  variability  of 
received  underwater  sound  can  improve  processing  and  handling  of  signals  of  interest,  for  example  the 
remediation  of  signal  degradation,  the  exploitation  of  intermittent  sonic  information,  or  the  exploitation 
of  favorable  propagation  conditions. 

OBJECTIVES 

The  main  goal  is  to  discover,  analyze  and  describe  mean  properties  of  acoustic  fields  and  time-varying 
acoustic  effects  caused  by  prominent  deterministic  propagation-altering  events  and  phenomena. 
Example  acoustic  effects  would  be  variation  of  horizontal  coherence  length,  variation  of  signal  level, 
variation  of  array  gain  (coherent  spatial  field  summation),  and  properties  of  the  (location-dependent 
and  time-dependent)  spatial  covariance  properties  of  acoustic  fields.  I  plan  to  study  specific  acoustic 
effects  caused  by  features  that  many  believe  are  the  most  influential:  water-mass  fronts,  internal-wave 
packets,  beam-like  internal  tides,  and  highly  structured  seafloors.  These  features  can  impart  non- 
Gaussian  form  to  the  acoustic  field  statistics.  The  seafloor  features  can  impart  spatially  non-Gaussian 
behavior,  with  example  features  being  sand  ridges,  slopes,  valleys  and  canyons. 

A  secondary  goal  is  to  more  fully  describe  the  propagation-altering  oceanic  phenomena  themselves. 
Working  towards  this  goal  is  closely  coupled  with  acoustics  research,  because  we  first  need  to  know 
which  features  are  the  more  relevant  for  acoustics,  and  to  what  detail  we  need  to  know  them.  This 
environmental  research  would  be  undertaken  with  acoustics  as  the  main  motivating  factor. 

APPROACH 

Although  a  lot  of  progress  has  been  made  in  describing  ocean  acoustic  variability,  numerous  factors 
motivate  continued  research,  among  them  the  high  complexity  of  the  ocean  environment,  the 
associated  complexity  of  the  propagated  sound,  and  the  lack  of  a  simple  conceptual  framework  for  the 
basic  properties  of  shallow-water  propagation  regimes.  The  research  direction  is  toward  such  a 
framework,  to  be  exploited  in  the  use  of  underwater  sound  in  shallow  water.  Our  knowledge  of 
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acoustic  field  patterns  in  shallow  water,  building  block  towards  a  useful  shallow-water  acoustic 
framework,  has  made  strong  advances  during  the  last  20  years  of  field  studies  and  theoretical  studies. 
We  wish  to  continue  this  progression,  establishing  the  relative  real-world  importance  of  various 
fluctuation-relevant  ocean  features,  how  the  features  interact,  and  how  the  sound  waves  move  through 
the  interacting  features.  In  this  approach,  models  of  the  environment  are  developed,  through  which 
sound  propagation  is  simulated  using  both  theoretical  and  purely  computational  methods,  including  our 
3D  parabolic  equations  solvers  (Lin  et  al.,  2013).  The  methodology  includes  building  environmental 
models  of  hierarchical  complexity,  from  which  the  salient  acoustic  processes  and  mechanisms  can  be 
identified  and  quantified. 

The  acoustic  simulation  tools  and  the  environmental  models  are  both  key  to  the  studies.  The  models  of 
the  variable  waveguide  features  can  be  both  idealized  and  realistic,  with  idealized  models  used  to  study 
parameter  dependence  and  uncover  basic  physics,  and  realistic  models  used  to  compare  computed 
and/or  theoretical  results  with  experiment.  Lying  between  these  two  activities  is  an  implicit 
comparison  of  results  from  the  idealized  and  realistic  models,  which  may  answer  the  question  “can 
insights  obtained  from  the  simple  models  be  used  to  understand  results  that  arise  from  the  realistic 


models”? 


The  acoustic  simulation  and  the  ocean  modeling,  and  the  analysis  of  output  will  be  run  by  the  PI  and 
by  Arthur  Newhall,  in  collaboration  with  other  members  of  the  WHOI  Ocean  Acoustics  and  Signals 
Laboratory.  At  this  time,  five  acoustic  studies  are  planned: 

(1)  Canyon  and  slope  acoustics:  Study  of  purely  geometrically  controlled  (no  ocean  dynamics) 

vs.  time-dependent  acoustic  field  features. 

(2)  Rough-seafloor  shelf  acoustics,  and  the  role  of  detailed  bathymetry  and  sub-bottom 

properties  in  the  presence  of  internal  waves. 

(3)  Comparison  of  realistic  vs.  idealized  ocean  model  inputs  to  acoustic  models. 

(4)  Evaluate  parameters  associated  with  array  signal  processing  and  adaptive  signal  processing. 

(5)  Canyon  experiment  planning:  (2015-16  for  an  anticipated  2017  or  beyond  experiment). 

WORK  COMPLETED 

Initial  calculations  have  been  examined  for  topic  numbers  one  and  four,  forming  a  parameter  study  for 
a  canyon  environment.  In  this  study,  horizontal  variability  in  the  acoustic  field  T(x,y,z)  is  quantified  in 
terms  of  the  horizontal  correlation  function  given  by 


where  the  averaging  indicated  by  the  brackets  can  be  computed  with  long  time  averages  for  stationary 
ergodic  processes,  and  may  include  spatially  averaging  over  x0  for  fields  with  homogeneous  statistics. 
This  statistic  was  examined  in  detail  by  Duda  et  al.  (2012)  for  internal  wave  effects  over  a  relatively 
flat  shallow  seabed,  and  we  move  here  to  an  entirely  different  regime.  Note  that  7^(0)  =  1.  The  point 
x  =  L  such  that  Rn(L)  =  exp(-l)  is  referred  to  as  the  horizontal  correlation  length  (or  coherence  length) 
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L.  The  spatially  lagged  product  in  the  numerator  tabulates  the  field  variability,  usually  dominated  by 
phase  differences  for  typical  ocean  acoustic  fields.  The  lagged  product  behavior  is  linked  to  array  gain 
for  additive  beamformer  processors,  which  use  sums  of  spatially  lagged  signals  rather  than  products. 
Gain  and  correlation  can  be  expressed  in  terms  of  one  another  for  certain  functional  forms  of  RN  (e.g. 

Cox,  1973). 

The  field  T(x,y,z)  to  be  inserted  into  the  formula  deserves  some  discussion.  For  a  fluctuating  process, 
the  time  mean  is  often  subtracted  from  each  point,  and  the  residual  forms  the  field  to  be  analyzed.  This 
would  provide  the  spatial  coherence  behavior  of  the  fluctuating  field,  but  the  mean  values  may  need  to 
be  known  and  utilized  to  make  the  best  use  of  an  extended  array,  and  may  not  be  available.  This  brings 
us  to  a  second  option  where  field-estimated  mean  values,  or  mean  values  fitting  an  a  priori  model  are 
used.  A  simple  a  priori  model  is  the  curved  or  spherical  wavefront  used  for  nearfield  beamforming.  A 
recent  paper  from  our  group  studied  apparent  acoustic  field  decorrelation  from  ocean  features  such  as 
front  intervening  between  a  sound  source  and  an  array  location,  which  would  impart  impart  a  steady 
but  unknown  structure  to  VF  (Lynch  et  al.,  2014).  If  known,  the  feature  and  the  structure  could  be 
accounted  for  to  achieve  gain  over  a  very  large  array  aperture  (i.e.  creating  a  very  large  effective  L ), 
but  currently  most  ocean  features  are  not  well  enough  mapped  to  be  accounted  for  in  this  way,  and 
actual  L  and  effective  L  are  often  relatively  small. 

After  completion  of  some  calculations  for  San  Diego-area  fluctuating  acoustic  fields,  we  are  now 
looking  into  the  implications  of  partial  array  coherence.  The  array  gain  AG  for  an  A-element  line  array, 
the  power  of  the  beamformed  signal  divided  by  the  power  of  the  beamformed  noise,  divided  the  signal 
to  noise  ratio  for  a  single  sensor,  is  10  logio  (N2 /N)  for  a  perfectly  coherent  signal  and  incoherent  noise. 
This  can  be  written  as  AG  =  SNRarray  /  SNRsensor,  where  SNRsensor  =  Ssensor  /  Nsensor  for  a 
single  sensor.  For  partial-array  coherence,  the  coherence  length  can  give  AG  for  an  arbitrarily  long 
array  by  substituting  N  +(L/Ay)  -  (L/Ay)  for  N  in  the  numerator,  where  Ay  is  the  inter-element  spacing 
in  meters,  so  that  AG=10  logio  ( 1  +\(L/ Ay)2  -  (L/ Ay)]/N).  The  quantity  L/Ay  is  the  number  of  effectively 
coherent  sensors  in  the  array.  The  quantity  AG  is  useful  because  it  can  be  used  to  estimate  SNRarray 
from  computations  of  L  made  from  3D  parabolic  equation  simulations.  In  log  units,  SNRarray  =  AG  + 
Ssensor  -  Nsensor.  SNRarray  is  tied  closely  to  the  detectability  of  a  particular  signal  of  interest.  As 
long  as  AG  +  Ssensor  exceeds  Nsensor,  which  is  the  ambient  noise  level,  the  system  is  usable.  Plots  of 
results  from  this  calculation  are  presented  in  the  Results  section. 

Additional  work  completed  in  these  first  months  consisted  of  revisions  to  two  manuscripts  to  appear  in 
the  IEEE  Journal  of  Oceanic  Engineering.  These  are  listed  at  the  end  of  this  report.  The  Lin  et  al.  paper 
describes  in  detail  the  measurement  and  modeling  of  sound  above  a  canyon  near  Taiwan,  with  the 
measurements  made  during  the  ONR-sponsored  Quantifying,  Predicting,  and  Exploiting  Uncertainty 
(QPE)  program  (Gawarkiewicz  et  al.,  2011;  Gawarkiewicz  and  Jan,  2013;  Newhall  et  al.,  2010).  The 
Emerson  et  al.  paper  describes  fluctuations  of  sound  energy  received  in  controlled  acoustic  studies 
performed  as  part  of  QPE,  and  implications  for  sonar  system  performance. 

RESULTS 

Figure  1  shows  results  for  a  3D  parabolic  equation  computational  example  of  array  performance 
prediction,  as  explained  above.  In  the  lowest  panel  AG  +  Ssensor  is  plotted,  in  dB,  for  a  700-Hz  sound 
source  radiating  down  Coronado  Canyon,  California.  The  AG  is  computed  from  L  estimates  (top 
panel)  made  using  many  scattered  imaginary  horizontal  test  arrays  of  100-m  meter  length,  with 
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element  spacing  Ay  of  3  meters.  Comparison  with  Ssensor,  which  is  plotted  in  the  center  panel,  shows 
that  the  array  provides  useful  gain  in  some  locations.  Figure  2  shows  results  for  a  similar  computation 
with  a  20-m  test  arrays.  Comparison  of  the  lower  panels  of  the  two  figures  show  that  the  different  size 
arrays  have  different  performance  characteristics,  with  processed  signal  level  exceeding  100  dB  in 
different  portions  of  the  plots,  with  the  short  array  performing  better  close  to  the  source.  In  the  very 
deepest  areas  studied  here  (upper  right  of  the  domain)  the  long  L  enhances  array  performance. 


60  m  depth  H  correlation  scale 


x  (km) 


Ssensor,  array  length  100  m  ;  set  1 


1  1.5  2  2.5  3  3.5 


x  (km) 


AG  +  Ssensor,  set  1 


1  1.5  2  2.5  3  3.5 

x  (km) 


■  110 

I100 

•90 


Figure  1.  3D  parabolic  equation  simulation  outputs  for  700-Hz  sound  going  down  Coronado 
Canyon  from  left  to  right  ( 150  dB  source  is  50  m  deep  at  fx,y]=[0,0J).  Black  lines  are  depth  contours 
withlOO-m  increment;  the  deepest  at  the  upper  right  is  700  m.  (upper)  The  horizontal  correlation 
length  L  in  they  direction,  calculated  at  60  m  depth.  The  calculation  is  made  with  a  100-m 
imaginary  aperture,  so  L  cannot  exceed  100  m.  (center)  The  average  signal  level  for  sensors  on  the 
imaginary  arrays,  (lower)  The  SNRarray  computed  as  described  in  the  text,  with  sensor  spacing  of  3 
m  assumed.  Comparison  of  the  lower  two  panels  shows  that  the  full  array  can  provide  useful  gain  in 
deep  water  at  the  upper  right,  but  provides  little  useful  gain  at  the  low  edge  of  the  domain,  which  is 

the  north  slope  of  Coronado  Canyon. 
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Figure  2.  The  same  signal  and  array  performance  predictions  as  shown  in  Figure  1  are  made  for  an 
imaginary  array  of  20-m  length  instead  of  100-m  length.  The  short  aperture  array  provides  gain 
near  the  source  ( compare  center  and  bottom  panels)  that  does  not  occur  for  the  longer  aperture. 

IMPACT/APPLICATIONS 

The  ultimate  performance  limitations  of  marine  acoustic  systems  are  set  by  the  properties  of  the  sound 
signals  that  propagate  within  the  fluctuating  three-dimensional  oceanic  environment.  Knowledge  of  the 
features  that  impact  sound  fields,  and  knowledge  of  how  the  features  affect  the  sound,  can  guide  design 
and  use  of  acoustic  systems. 

RELATED  PROJECTS 

The  PI  also  has  a  MURI  grant  managed  by  ONR  that  is  closely  related.  That  project  is  working  toward 
coupled  modeling  of  the  ocean  environmental  and  acoustical  modeling.  Information  can  be  found  at 
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http://www.whoi.edu/sites/IODA.  The  PI  also  has  a  National  Science  Foundation  Physical 

Oceanography  Program  grant,  with  Co-PIs  Y.-.T.  Lin  and  Pierre  Lermusiaux,  to  study  ocean  internal 

tide  dynamics. 
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